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Abstract. 

Employing a local moment approach to the periodic Anderson model within the 
framework of dynamical mean-field theory, direct comparison is made between theory 
and experiment for the d.c. transport and optical conductivities of paramagnetic heavy 
fermion and intermediate valence metals. Four materials, exhibiting a diverse range of 
behaviour in their transport/optics, are analysed in detail: CeB^, YbAl^, CeAl^ and 
CeCoIn^. Good agreement between theory and experiment is in general found, even 
quantitatively, and a mutually consistent picture of transport and optics results. 
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1. Introduction. 

Heavy electron materials have long been the subject of extensive investigation, for 
reviews see e.g. [1-7] . Yet in many respects even their 'normal' paramagnetic phase, be it 
metallic or insulating, has eluded a unified microscopic description on all experimentally 
relevant temperature (T) and/or frequency (u) scales. The canonical theoretical model 
here is of course the periodic Anderson model (PAM). Within the general framework of 
dynamical mean- field theory [8-11] we have developed in the preceding paper [12] (here 
referred to as I) a non-perturbative local moment approach to paramagnetic metallic 
phases of the PAM, with a focus on d.c. transport and optics; following earlier work 
on T = dynamics [13] and on Kondo insulators [14, 15]. The primary emphasis of I 
is the Kondo lattice regime relevant to strong correlated heavy fermion (HF) metals. 
Dynamics/transport on all relevant (u, T) -scales are encompassed, from the low-energy 
behaviour characteristic of the lattice coherent Fermi liquid, through incoherent effective 
single-impurity physics to non-universal high-energy scales. The underlying theory 
is not however restricted to the Kondo lattice regime, enabling it also to handle e.g. 
intermediate valence (IV) behaviour. 
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The present paper is an attempt to provide at least a partial answer to the question: 
to what extent are the optical and d.c. transport properties of HF and related materials 
captured by the PAM and our theory for it? That clearly requires direct, quantitative 
comparison between theory and experiment, which is our purpose here. Specifically, we 
consider in detail the transport and optics of four materials: CeB 6 , YbAl 3 , CeAl 3 and 
CeCoIn 5 , all HF metals with the exception of the IV compound YbAl 3 , and exhibiting 
a diverse range of behaviour in their transport and optical behaviour. The materials 
are analysed on a case by case basis in §s 3-6, following a discussion (§2) of relevant 
issues involved in making the comparison; and we believe it fair to claim that the theory 
provides a striking account of experiment. 

2. Background issues 

The Hamiltonian for the PAM is given by equation (2.1) of I, and its physical content is 
simple: a single correlated /-level in each unit cell hybridizes locally to an uncorrelated 
conduction band. The model is moreover specified by only four 'bare '/material 
parameters rendering it minimalist in terms of comparison to experiment — the more 
so when one recalls that it encompasses regimes of behaviour as diverse as strongly 
correlated heavy fermion metals and Kondo insulators, intermediate valence, and weak 
coupling. The dimensionless bare parameters are U, V, e c and tf (in units of the 
conduction electron hopping, t* = 1), with U denoting the local /-level interaction 
strength and V the local one-electron hybridization coupling an /-level to the conduction 
band. The energy of the local conduction orbital, e c , determines the centre of gravity of 
the free conduction band relative to the Fermi level (and thereby the conduction band 
filling, n c ); and e/ denotes the /-orbital energy. An equivalent parameter set is U, V, 
e c and rj, where rj = 1 + 2e/ /U specifies the /-level asymmetry. 

The non-interacting limit of the model (U = 0) is certainly trivial. But in 
that case — for all T — the d.c. resistivity of the metallic state vanishes, and the 
optical conductivity contains no absorption below the direct gap save for a 5-function 
Drude contribution at u — 0, see I. That this behaviour bears scant comparison to 
experiment reflects the fact that the essential physics is driven by scattering due to 
electron interactions. It is of course the latter, and the resultant many-body nature of 
the problem, that renders the PAM non-trivial. 

In considering d.c. transport, the first requirement in comparing theory to 
experiment is thus to extract the contribution (jt)^(T)) to the measured resistivity 
(p(T)) that isolates the interaction contributions from those of phonons and static 
impurity scattering. This is given, ideally, by 

«CO = {p(T) - p(0)} - p ph (T). 

The first term removes the residual (T = 0) resistivity, and hence impurity scattering 
contribution on the assumption that the latter is T-independent. The second removes 
the contribution from phonons; usually taken in practice (as assumed in the following) 
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to be the resistivity of the non-magnetic homologue compound with the magnetic ion 
Ce (or Yb) replaced by La (or Lu), on the assumption that interactions in the latter 
are negligible. For most systems the phonon contribution is generally negligible for 
T < 50K or so. 

While this prescription is straightforward in principle, we first deal with a 
complication that can arise in practice. Experimentally it is the resistance that is 
measured directly. To convert to a sample independent resistivity requires spatial 
dimensions to be known with reasonable precision. That does not pose a problem with 
large crystals, but may do for small samples. Examples arise in the literature where 
reported p(T)'s from different groups differ quite significantly; we encounter one such 
in the case of CeCoIn 5 considered in §6. Absolute resistivities would thus be related 
to measured values by e.g. a'p(T) and a"p p h(T), where a' and a" denote experimental 
'mismatch' factors. For comparison to theory (equation (2.2) below) we require however 
only the relative factors a' /a", to which end we replace the above by 



of I for the PAM (where p mag was denoted simply by p). The system is generically 
characterized by a low-energy coherence scale u>l = ZV 2 , with Z the quasiparticle 
weight or inverse mass renormalization factor. This scale is a complicated function of 
the underlying bare parameters, see e.g. [13] and refs therein. But in the strong coupling 
Kondo lattice regime ojl is exponentially small (because Z is). In consequence, p mag (T) 
exhibits scaling in terms of u>l, i-e. is of form 



with the temperature dependence encoded in T/cul, independent of the interaction 
strength U and hybridization V (a is a trivial overall scale factor, a = 1/gq in the 
notation of I); the scaling holding for (any) fixed e c and rj. The scaling resistivity is 
moreover only weakly dependent on 77, and for T/oul > 1 — 5 or so is in fact independent 
of e c , reflecting the crossover to incoherent effective single impurity physics (as detailed 
in I, see e.g. figure 8). This enables direct connection to experiment, via the extent 
to which the scaling form equation (2.2) captures the T-dependence of experimental 
resistivities; and indeed also their pressure dependence, for although ojl will change 
with pressure the scaling behaviour should remain intact (we consider this explicitly in 
the case of CeAl 3 , §5). Success in this regard also enables lo l to be determined directly, 
there being no hope of calculating it 'ab initio' with any meaningful accuracy. Finally, 
we add that comparison of theory/experiment proceeds along the same lines away from 
the asymptotic Kondo lattice regime, in dealing e.g. with intermediate valence materials 
or intermediate/ weak coupling compounds. Here a full parameter set e c , rj, U and V 
must in general be specified; but p ma g(^) can always be cast in the form equation (2.2), 
and appropriate comparison can be made. 





(2.2) 
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Although our comments above focus on static transport, a central purpose of the 
paper is also to make direct comparison between experiment and theory for optics, on 
all experimentally relevant frequency and temperature scales. Given prior analysis of 
d.c. transport, the underlying model/material parameters are known, either wholly or 
in part. There is then little room for manoeuvre; the resultant theory either captures 
the optical behaviour or not, providing a further and quite stringent test of the model's 
material applicability. 

The main physical effect omitted in the PAM itself is that of crystal electric fields 
(CEFs). The atomic levels of e.g. Ce 3+ , 2 Fs, are generically split into three doublets, 
the excited levels lying above the ground state by Ai, A 2 . For sufficiently low T only 
the lowest level matters, it being this alone the PAM seeks to capture. Although 
the material specific A« are usually larger than uol (itself typically ~ 10 — lOOif), 
they often lie in the interval < 300/T. Their qualitative influence on d.c. transport is 
clear, for additional conduction channels are opened up in accessing higher CEF levels 
with increasing T. Quantitatively however, the effect is hard to gauge a priori, its 
magnitude naturally depending on how effectively the higher CEF levels couple to the 
conduction band. Where present and effective, we can expect to see it as a decrease in 
the experimental resistivity below that predicted by the 1-channel PAM, the onset of 
the deviation appearing at T ~ A^ We add here that while the role of CEF effects 
has been studied in the context of smg/e-impurity Anderson/Kondo models [16-20], 
with application to lattice-based systems for sufficiently high temperatures where lattice 
coherence can be neglected, we are not aware of corresponding work in the context of 
lattice-fermion models. 

Finally, we note that our comparison of theory/experiment for Kondo insulators 
in [15] was free from most of the considerations above. There, theory was compared 
directly to the experimental p(T), for several reasons. First, in contrast to the case 
of heavy fermion metals, resistivities of the non-magnetic homologues are sufficiently 
small compared to those of the Kondo insulating material that p p h(T) in equation (2.1) 
can be neglected with impunity. That in turn means that any sample geometry factor 
a 7^ 1 in equation (2.1) can simply be absorbed into the overall scale factor a (equation 
(2.2)), so any lack of precision in obtaining p(T) is immaterial. The role of impurities is 
also different in Kondo insulators. In metals this arises from static impurity scattering, 
presumed to be T-independent and generating the finite residual resistivity p(0); which 
is thus subtracted out as in equation (2.1). In Kondo insulators by contrast p(T) diverges 
as T — > 0, reflecting the insulating ground state. This occurs even in the presence of 
localized impurity states in the T = insulating gap; which generate conduction by 
variable-range hopping, generally operative over a narrow T-interval (e.g. < 8K for 
SttiBq [15,21]), and whose net effect on p(T) at temperatures above this interval (where 
direct comparison is made to experiment) is usually sufficiently small to be neglected. 

We turn now to comparison with experiment for the metallic heavy fermion and 
intermediate valence materials considered. 
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3. CeB 6 

Among the rare-earth hexaborides, the Kondo insulator SmB^ and the heavy fermion 
compound CeB 6 have been investigated for many years [1, 7, 22-26]. The former 
was considered in our recent work [15] on the particle-hole symmetric PAM. There 
it was shown that a single low-energy (indirect gap) scale underlies the temperature 
dependence of both the static and optical conductivity, and (with minimal input of bare 
material parameters) the frequency dependence of the optics as well. Here we consider 
its metallic counterpart CeB 6) likewise a cubic system [22]. At the lowest temperatures, 
various antiferromagnetic phase transitions occur between T ~ 1.6 — 3.3/T [22,24,26]; 
above (and below) which the paramagnetic phase arises, 'Phase-F for T > 3.3i^. 

The relative ease with which large, clean single crystals of CeB e can be grown (see 
e.g. [26]) has been a motivating factor in its investigation. The d.c. resistivity has been 
measured by several groups (e.g. [22, 26]) and since the sample quality is in general 
good, the residual resistivity is very small. In addition, the relatively large crystal sizes 
imply small errors in sample geometry, hence the factor a in equations (2.1) can safely 
be assumed to be 1. The phonon contribution to the resistivity is as usual taken as 
the resistivity of the non- magnetic homologue LciBq [22], the latter having the same 
lattice structure with similar lattice parameters and phonon dispersion as CeB e (from 
inelastic neutron scattering [22]). The experimental magnetic resistivity p^i g {T) is then 
readily obtained [22], and shown in figure 1 (we add that the phonon contribution 
kicks in only for T > 50K, so the magnetic resistivity essentially coincides with the 
raw resistivity at lower temperatures). A small kink arises in p^fgCO at T ~ 3.3K, 
reflecting the transition from paramagnetic phase I to phase II [22,26] mentioned above. 
The experimental resistivity is seen to exhiCoverletter fig4.eps iopartlO.clo paper.tex 
ReplyLetter bit the classic 'shape' for a strongly correlated HF metal, increasing from 
zero at T = and going through a coherence peak at T ~ 4K, before decreasing through 
a small log-linear regime (similar to that shown in figure 6(a) of I) to a shallow minimum 
at T ~ 375K; increasing thereafter at higher temperatures where it shows conventional 
metallic behaviour. 

To compare p^ g (T) to theory, first recall from I that the asymptotic scaling 
resistivity Pmag(^) is a universal function of T/ul for fixed e c and r\ (being independent 
of the local interaction U and hybridization V); see e.g. figure 8 of I. For a chosen 
(e c ,v) the scaling p n mg(T) is then straightforwardly superposed onto p^l g {T) with an 
appropriate rescaling of the temperature and resistivity axes, thus enabling the low- 
energy scale ujl to be determined. The asymptotic scaling p mag (T) is shown in figure 
1 (dotted line) for a moderate e c = 0.5 and rj = 0, with the resultant coherence scale 
thereby found to be ojl — 5.5/T. It is seen to capture the experimental resistivity up 
to T ~ lOOfr, but above this it deviates below experiment, continuing as it must to 
decrease monotonically (see figure 8 of I) and hence lacking the minimum occurring 
experimentally for CeB e at T min ~ 375K, i.e. T min = T m - m /uj L ~ 70. 

As discussed in I (figure 8 inset), the latter behaviour is physically natural and 
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readily encompassed theoretically. No real HF material is in the universal scaling regime 
'for ever' — it must be exited sooner or later with increasing T. Deviation of p mag (T) 
from its asymptotic scaling form at sufficiently high temperatures signifies the onset of 
non- universality; and the location of the non- universal minimum in p mag (T) provides 
an opportunity to identify the ratio U/V 2 of effective bare material parameters, which 
will be helpful in making a prediction for the tu-dependence of the optical conductivity 
a(u; T). Specifically, for the chosen e c = 0.5, 77 = 0, we find that the theoretical p ma g(^) 
indeed has a minimum at T min ~ 70 (as in experiment) for U /V 2 ~ 12. This is illustrated 
in figure 1, where we show the resultant theoretical p mag (T) for U = 4.75, V 2 = 0.4 (solid 
curve) and U = 2.45, V 2 = 0.2 (dashed curve). The two p m ag(7 1 )'s barely differ from the 
asymptotic scaling resistivity for T < 100K, are essentially coincident with each other 
across the entire T-range, and each possesses a weak minimum at T ~ 375i^. Except 
naturally for a small neighbourhood around the low-temperature phase transitions 
(which the theory does not seek to address), the resultant theoretical p ma g(^) is seen 
to be in rather good agreement with experiment, as evident further from the inset to 
figure 1 where the corresponding d.c. conductivity a mag (T) is shown. A CEF excitation 
is known to occur at ~ 2.5meV (30K) [27] (with a second lying at a much higher energy, 
46meV [27]); but as judged from the above comparison this appears to play a minor 
role in the d.c. transport itself. 



100^ 




Figure 1. Comparison of experimental p^i g (T) (circles) for CeB^ [22] to theory, on 
a log-linear scale. The solid curve shows /O mag (T) for U = 4.75, V 2 = 0.4, while the 
dashed curve is for U = 2.45, V 2 = 0.2; both theory sets have common e c = 0.5, rj = 
and the same 00 l — 5.5K. The dotted curve shows the asymptotic scaling resistivity. 
Inset: the d.c. conductivity cr mag (T) = l/p mag (T) on a log-log scale. 



To determine the optical conductivity on all a>scales, including non-universal 
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frequencies, requires U and V 2 to be separately specified as discussed in I. That is 
clearly not provided by the above analysis, but a reasoned prediction for a(u;T) 
may nonetheless be made. For fixed U/V 2 ~ 12 as above, we find that varying U 
across the range ~ 2.5 — 5 produces only a modest change in both the uj- and In- 
dependence of a(u;T), leading in particular to a direct gap absorption lying in the 
interval ~ 200 — 300 cm -1 ; as well as a quasiparticle weight Z on the order of 10" 2 
that is consistent with the effective mass (m*/m e =) m* ~ 100 deduced from the 
specific heat coefficient 7 = 250 mJmol^K^ 2 measured for Phase-I of CeB 6 [22,26]. 
In otherwords, on the assumption that the reasonably wide [/-range above encompasses 
the behaviour of CeB$, the optics are relatively insensitive to the precise value of U. 
Coverletter fig4.eps iopartlO.clo paper.tex ReplyLetter 

To that end we show in figure 2 the predicted optical conductivity for U = 2.45, 
V 2 = 0.2, with the u;/u;L-dependence of the theoretical cr(u; T) converted to uj in cm" 1 
using uj l = 5.5K(= 3.8cm _1 ) deduced from the above analysis of the d.c. resistivity; 
and with cr(u; T) shown for a range of temperatures from 1.1K to 660K. As uj — > 0, the 
T-dependence of the dynamical conductivity follows the d.c. values shown in the inset to 
figure 1. At the lowest temperature T — 1.1K ~ ul/5 an emergent low- frequency Drude 
absorption on frequency scales uj < 1 cm" 1 is evident in cr(u; T) (see also figure 9 of I), 
separated by a clear optical pseudogap centred on uj ~ 10 cm" 1 from a strong direct 
gap absorption centred on uj ~ 250 cm" 1 . The Drude absorption is rapidly suppressed 
on increasing T, being all but dead by T ~ ujl = 5.5K; while the optical pseudogap is 
progressively 'filled in' on temperature scales set by ujl- The direct gap absorption is 
largely unaffected by temperature until T > (5 — 10)uj l or so; but is significantly eroded 
by T ~ HOfT and in essence destroyed by room temperature. 

The optical conductivity of CeB Q has in fact been measured by Kimura et al [23], 
but only in the frequency interval 50meV ~ 400cm" 1 to lOeV ~8x 10 4 cm _1 , and at a 
temperature of 300fT. No strong direct gap absorption was observed, itself suggesting 
that the absorption occurs at frequencies on the order of uj ~ 200 cm" 1 ~ 300/T or 
less; instead a broad, featureless and monotonically decreasing spectral lineshape was 
found. This is of course consistent with the prediction from figure 2 above that the direct 
gap peak is almost completely washed away at room temperature. In order to observe 
non-trivial uj- and T-dependence in the optical conductivity of CeB Q , we thus suggest 
the frequency domain be extended down to ~ 10 cm" 1 (or lower to observe the Drude 
absorption), and that experiments be performed at considerably lower temperatures 
such as those shown in figure 2. 

4. YbAh 

Our main focus in I and [13] has been the strongly correlated heavy fermion regime 
where the /-level e/ < lies well below the Fermi level, with e/ + U 3> well above it, 
such that the /-electrons are essentially localized, rif — > 1. The underlying local moment 
approach is not however restricted to the Kondo lattice regime, and in particular can 
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Figure 2. Predicted optical conductivity of CeB§: <j(gj;T) vs. uj in cm 1 at 
temperatures T = 1.1 A - , 5.5K, 33if, 55^, llOif, 330-fT, 550X and 660K. 



also readily handle intermediate valence (IV) behaviour. In this case, depletion for 
example of nj from unity reflects the fact that ej lies relatively close to the Fermi level, 
such that the appropriate r\ — 1 + 2e//C/ regime is rj ~ 1 . 

The compound YbAl 3 , which crystallizes in a simple cubic Cu 3 Au-tjpe structure 
and does not order magnetically [28] , provides a prime example of IV behaviour. Figure 
3 shows the experimental d.c. resistivity of YbAl 3 , and its non-magnetic homologue 
LuAl 3 , at ambient pressure (data from [28] with the tiny residual resistivity subtracted). 
The resultant magnetic resistivity p^ g (T) is also shown in figure 3, obtained simply by 
subtracting the resistivity of LuAl 3 from that of YbAl 3 (the samples [28] are high quality 
single crystals, hence the factor a in equation (2.1) is taken as unity). p^ g (T) is seen to 
increase monotonically with T, lacking the coherence peak seen in strongly correlated HF 
materials. This behaviour is characteristic of IV [29], as too is e.g. the low/moderate 
effective mass in the range m* ~ 15 — 30 inferred from dHvA [30], optical [31] and 
specific heat [32] measurements; and the Yb mean valence z v = 2 + rif is estimated 
experimentally as z v ~ 2.65 — 2.8 [32,33]. 

In the strongly correlated HF regime, physical properties such as p mag (T) exhibit 
scaling as a function of T/ul, independent of the interaction U and hybridization V as 
detailed in I; occurring formally for all T/ujl in the asymptotic strong coupling limit, 
and in material practice over a significant albeit naturally finite T/ujl range as seen 
above for CeB§ (and in [15] for Kondo insulators). This is not by contrast the case 
in the IV regime, and neither is it to be expected. Here the 'full set' of bare/material 
parameters 77, e c , U and V must be specified. To compare theory to experiment for YbAl 3 
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Figure 3. The d.c. resistivity of YbAlz (diamonds) and LuAlz (point-dash line), 
from [28] with the small residual resistivity subtracted. The experimental p^ g (T) 
(circles) is obtained by subtracting p(T) — p(0) for L11AI3 from that for YbAl^. The 
theoretical p ma , g {T) (dashed line) is obtained for e c = 0.5, rj = 1.2, U = 4.9, V 2 = 0.8, 
and superimposed on p^[ g (T) with lol — 254K. The level of agreement between theory 
and experiment is clear. Inset: the low temperature behaviour; including (dotted line) 
a fit to the asymptotic Fermi liquid form p mag (T) oc T 2 , which is seen to persist in 
both experiment and theory up to T ~ 30K . 



we consider a moderate e c = 0.5 with rj = 1.2, U = 4.9 and V 2 = 0.8; corresponding 
to a modest interaction strength U/irA = 1.4 (with A = 7rV 2 po(— e c ) as in I) and an 
6f — 0.49 close to the Fermi level (e//Ao = 0.44). The resultant /-band filling is found 
to be nf ~ 0.65 (z v ~ 2.65), consistent with the mixed valence nature of YbAl 3 ; and 
a quasiparticle weight Z ~ 0.05 is found, implying an effective mass m* ~ 20 that is 
likewise consistent with experiment as above. 

The theoretical d.c. resistivity p ma , g (T) vs. T/ujl for these parameters has been 
determined, and superposed onto the experimental p^ g (T) in the usual way. The 
resultant low-energy scale is found thereby to be = 254i\~, and comparison between 
theory and experiment is shown in figure 3. The agreement is clearly excellent, for 
all temperatures. The inset to figure 3 shows the low temperature behaviour on an 
expanded scale, together with a fit (dotted line) to the T — > Fermi liquid form 
Pmag(r) oc T 2 . As known from experiment [32], and seen also in the present theory, 
this asymptotic form is seen to persist up to a temperature £ T"fl'~ 30K that is an order 
of magnitude lower than ojl ~ 254K; and about which fact we make two brief comments. 
First, to emphasise that this arises naturally here within a theoretical approach to the 
periodic Anderson model itself (with no appeal e.g. to calculations based upon a single- 
impurity Anderson model [32]). Second, the quantitative distinction between 'Tfl an d 
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T [K] 

Figure 4. Theoretical p(T) for YbAl^ (dashed line) up to 800^, compared to the 
experimental results of [28] (diamonds) up to T = 300K and of [34] from 300 — 750K 
(triangles), as detailed in text. 

u>l is in our view to be anticipated; for uj l is the natural low-energy scale in the problem, 
and since p(T) oc T 2 is the asymptotic T — > behaviour, we would as such expect it to 
arise only for T <C uo L . 

While the experimental data shown in figure 3 extend up to T ~ 300i^, the 
theoretical PmagiT) may be used to predict the full p(T) for YbAl 3 over a much larger 
temperature interval. To that end we simply calculate p ma ,g(T) at higher T, and add 
to it the resistivity of LuAl% representing the phonon contribution, itself extended to 
higher T by linear extrapolation of the Lu data in figure 3 (point-dash line, and which 
extrapolation is clearly warranted). The resultant p(T) is shown in figure 4 out to 
T = 800/T; and we note that it continues to increase monotonically above 300i^, 
precluding as such the occurrence of a coherence maximum at a T in excess of that 
shown in figure 3. Experimental results for p(T) up to T ~ 750/T have in fact been 
reported [34]. This data does not appear to be quite as clean as that of [28] in the 
interval 50/T < T < 300iT, but for T < 50fT the p(T) from [34] collapses very well onto 
that of [28] (considered above) with an overall p-axis rescaling factor a = 1.2 (equation 
(2.1)). Taking this a, the resultant p(T) from [34] is shown in figure 4 in the temperature 
interval 300 — 750K; and is seen to agree well with the theoretical result. 

We turn now to the optical conductivity of YbAl 3 , which has only recently been 
measured at infrared frequencies and below [31]. This is reproduced in the top panel 
of figure 5, from which three key features are evident [31]: the low-frequency Drude 
response characteristic of the free carriers; a depleted pseudogap occurring at low 
temperatures (at u ~ 20 meV) and flanked on its right by a shoulder at uj ~ 50—60 meV; 
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followed by the broad, strong direct gap (or mid-infrared, mIR) peak centred near 
uj ~ 250meV. 
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Figure 5. Top panel: experimental optical conductivity of YbAl% [31]. Bottom panel: 
theoretical cr(w; T) for the same parameters as in figure 3; at the same temperatures 
and on the same frequency scale as in experiment. The vertical arrow in both panels 
indicates the position of the theoretically predicted shoulder, at u) ~ 2ll>l ~ 50 meV. 
The insets show the thermal evolution of er(aj; T) for higher temperatures. The overall 
experimental lincshapc including the pseudogap, shoulder and mIR peak, as well as 
their thermal evolution, is well reproduced by theory. 



The bottom panel in figure 5 shows the theoretical optical conductivity at the same 
temperatures and on the same frequency scale as experiment, obtained using the same 
bare parameters employed in figures 3,4, and with u L = 254K as deduced above; in 
otherwords with no additional input other than that inferred from the d.c. transport 
comparison. Barring a mismatch in the relative intensity of the direct gap peak on 
which we comment below, the agreement between theory and experiment is seen to be 
rather good in terms of the pseudogap structure, the shoulder, the position and width 
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of the mIR peak, and the thermal evolution of the optical conductivity. The shoulder in 
particular merits comment, being a distinct spectral feature that is thermally destroyed 
with increasing temperature [31]. It has been speculated in [31] (not unreasonably) 
that its origin may lie outside the scope of the PAM and/or be material specific. The 
present results however suggest to the contrary. Theoretically, we find that the existence 
of the shoulder shown in figure 5 is not specific to the particular set of bare material 
parameters considered, but rather characteristic of IV in general terms — a clear optical 
shoulder arising in the vicinity of u ~ 2 — 3u>l as the bare parameters are varied over 
quite a significant range. Its origins reflect the underlying behaviour of the /-electron 
self-energy, and we do not therefore have a simple physical explanation for it; but neither 
do we doubt its generic occurrence. For the case of YbAl 3 , with uj l ~ 250-fT deduced 
from d.c. transport as above, the shoulder thus lies at u ~ 50 meV as indicated by the 
vertical arrows in figure 5 and agreeing rather well with experiment. 

Finally, we comment on the mismatch in the vertical (intensity) axes in figure 5. 
The theoretical a(u;T) represents of course the conductivity in the absence of phonons, 
which inevitably introduces a certain mismatch between theory and experiment as 
uj — > at temperatures high enough for the phonons to kick in; but the phonon 
background would not extend up to mIR frequencies, so the intensity of the direct 
gap peak in theory and experiment should agree. The fact that it does not could 
obviously mean that a different bare parameter set may be more appropriate than that 
identified here. Alternatively, the issue may be experimental. The cr(iv; T) spectra are 
obtained [31] via Kramers-Kronig transformation of the total reflectance between 7meV 
and 35eV, with a Hagen-Rubens formula used for low-energy extrapolation; and with 
seCoverletter fig4.eps iopartlO.clo paper.tex ReplyLetter veral spectrometer sources 
employed in different frequency intervals, which require matching using appropriate 
constant factors (see e.g. [35]). This is intricate, and there is undoubtedly the possibility 
of error in determining the absolute intensity of the direct gap. A surface impedance 
probe such as the one used in [35] might be able to resolve the matter; while if the issue 
is not experimental then a more extensive scan of the bare parameter space is required. 
That notwithstanding, however, the present theory does appear to provide a remarkably 
consistent description of both transport and optics in YbAl 3 . 

5. CeAh 

The classic system CeAl% has long been subject to extensive investigation, see e.g. 
[1,35-42]. In contrast to its Yb cousin considered above, CeAl 3 is a prototypical 
heavy fermion material; as attested for example by the large specific heat coefficient 
7 ~ 1.4 x 10 3 mJmol~ l K~ 2 [39] and corresponding effective mass m* ~ 700 [35,39]. 

A helpful starting point for comparison of theory to experiment is a rough knowledge 
of the parameter regime to which the system belongs. In the case of CeAl 3 we 
can glean such information from the unusual behaviour of the experimental optical 
conductivity [35], which is shown in the top panel of figure 7 below. Typical features 
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found in optical lineshapes of HF systems are the Drude peak at low frequencies, followed 
by a pseudogap, with a strong direct gap (or mIR) absorption at higher frequencies (see 
e.g. figures 2,5,10). In the case of CeAl^ however, while the Drude peak is clearly 
evident at low T, a distinct pseudogap and the mIR peak are absent. In fact there is 
little frequency dependence beyond the Drude absorption range. 

From our theoretical work in I we know that increasing e c , and hence reducing 
the conduction band filling n c , acts to diminish the pseudogap (see inset to figure 10 
of I); relatedly, it also tends to suppress the direct gap absorption, suggesting that a 
suitably large e c is required for CeAl 3 . However this by itself is not sufficient: a large 
e c , but with a modest hybridization V between the /-levels and the conduction band, 
can still give rise to a distinct mIR absorption. A suitably large hybridization thus 
also seems necessary to suppress strongly the direct gap absorption. Unusually large 
CEF parameters found in inelastic neutron scattering experiments have in fact also been 
attributed to a large hybridization V [40], supporting this inference. And from the large 
effective mass mentioned above, we know the system is in the strong correlated Kondo 
lattice regime, requiring a significant interaction strength U . 

The picture of CeAl 3 thus suggested is of a system with low conduction band filling 
n c (large e c ), a significant hybridization V, and strong local interactions. We have 
investigated this regime in some detail, and for comparison to transport and optical 
experiments on CeAl^ will consider explicitly the following material parameters: e c = 1.5 
and i] = (results are quite insensitive to r) in the Kondo lattice regime), together with 
V 2 = 1.4 and U = 6.9 (or U/ttAq ~ 8.4 ^> 1 signifying strong correlations, where 
A = nV 2 p (— e c ) as in I). Additional specification of U and V 2 is of course required in 
order to consider the optical conductivity on all frequency scales, whereas e c and rj alone 
suffice to determine the d.c. resistivity in the scaling regime. With these parameters we 
find a quasiparticle weight Z ~ 1.6 x 10~ 3 and hence an effective mass m* ~ 625, in 
good agreement with m* ~ 690 deduced from specific heat measurements [35,39]. The 
resultant conduction band filling is n c = 0.17; while the /-level occupancy rif = 1.0, 
consistent with the firmly trivalent (z v = 4 — nf) nature of the Ce ion expected in the 
HF regime. 

We first consider d.c. transport measurements. These have been obtained by several 
groups [35,37-39,41], which we find in general concur well on subtraction of appropriate 
residual resistivities (so a = 1 is taken in equation (2.1)). We choose to compare 
explicitly to the p(T) data of [37], which is shown in figure 6 (squares) at ambient 
pressure with the residual resistivity subtracted out; together with the corresponding 
Pmig(T) (circles) obtained by further subtracting the resistivity of LaAl 3 . Figure 6 
also shows corresponding results for p^ g (T) obtained at a pressure P = 0.4 GPa [38] 
(triangles), using the same sample as in [37]. The theoretical p mSug (T) is calculated 
and scaled onto the ambient pressure p^i g {T) in the usual way, leading thereby to a 
coherence scale of Ul ~ 33K. From figure 6, comparison between theory and experiment 
is seen to be good up to T ~ A5K, beyond which the experimental p e ^i g {T) drops much 
more rapidly with further increasing T. This is natural, for inelastic neutron scattering 
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Figure 6. d.c. resistivity of CeAl^ at ambient pressure (squares, from [37] with the 
residual resistivity subtracted); and the corresponding magnetic resistivity PmJgCO 
(circles) [37]. The theoretical p m ag(T) for e c = 1.5, y = is shown by the 
solid/dotted line, the resultant coherence scale being uj l — 33if. The agreement 
between theory and experiment is good for T < 45JT, beyond which deviations 
naturally occur due to the presence of two crystal field split levels at 6.1meV (71K) 
and QAmeV (75K) [40], marked by arrows. Triangles denote p2 g (T) for CeAlj, at a 
higher pressure P = 0.4 GPa [38]. The dashed line shows the same theoretical curve 
as used for comparison to the ambient pressure data, but with simple rescaling of the 
axes, leading to w_l(0.4 GPa) = 71K. 

experiments [40] show two higher crystal field levels occuring at almost the same energy, 
QAmeV (70. 8K) and QAmeV (7A.2K) above the ground state. Marked in figure 6, these 
are accessed thermally as T approaches ~ 70K, and provided they couple effectively 
to the conduction band (as appears to be the case here) the two additional conduction 
channels diminish the resistivity significantly. 

We now comment briefly on the high pressure magnetic resistivity shown in figure 6. 
The coherence scale ujl = ZV 2 itself naturally varies with pressure, but in the strongly 
correlated Kondo lattice regime the T-dependence of the magnetic resistivity should 
depend universally on T/uol alone (§2). The same theoretical p ma ,g{T) employed for 
comparison to the ambient pressure data should thus, with mere rescaling of the axes, 
account for p^ g (T) at P = 0.4 GPa (again up to the temperature at which the CEF 
excitations kick in). That indeed it does is shown in figure 6 (dashed line); the resultant 
coherence scale being found to be tu L (P = 0.4 GPa) ~ 71K. 
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Figure 7. Top panel: experimental optical conductivity of CeAls at the temperatures 
indicated, from [35]. Bottom panel: theoretical a(co;T). Discussion in text. The 
vertical arrow denotes the location of the direct gap, where a mIR peak is usually 
found for other HF systems. 



Finally we turn to the optical conductivity, the experimental results [35] shown in 
figure 7 (top panel) being obtained from reflectance spectroscopy (lines) and surface 
impedance measurements (points). The bottom panel shows the theoretical a(w,T) 
obtained with the bare parameters specified above, converted to cm -1 using ujl = 33K 
inferred above from static transport. The theory evidently captures the unusual optical 
characteristics of CeAl 3 , and agreement with experiment is seen to be rather good. A 
strong Drude peak, and a very shallow pseudogap, are seen at the lowest temperatures 
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in theory and experiment; the Drude absorption collapsing on a temperature scale of 
~ 10K, above which very little T dependence is found across the entire frequency range. 
Most significantly, no distinct direct gap/mIR peak arises, the spectrum being largely 
featureless (the small feature visible at ~ 10 4 cm _1 in the theoretical o~(uj; T) occurs on 
the effective bandwidth scale, and its intensity diminishes further with increasing U). 
The nominal location of the direct gap itself can however be determined theoretically 
frCoverletter fig4.eps iopartlO.clo paper.tex ReplyLetter om the renormalized band 
structure underlying the present theory (as considered in I, see figure 11). It is found 
to lie at uj ~ 10 3 cm _1 as marked in figure 7, albeit that no sharp absorption occurs in 
its vicinity. 

6. CeCoIn 5 

This recently discovered [43], moderately heavy fermion compound crystallizes in a 
tetragonal structure consisting of alternate layers of Celn 3 and CoIn 2 . It superconducts 
below T c ~ 2.3K (the highest transition temperature of all known HF systems at 
ambient pressure [43]), and is paramagnetic for T > T c . The experimental resistivities 
p(T) determined by four different groups [43-46] are shown in figure 8. That they differ 
widely presumably reflects intrinsic difficulties in measuring the dimensions of relatively 
small samples. Their basic equivalence is however seen by taking the results of one 
as a reference and rescaling the y-axis for each of the remaining data sets. With this, 
as shown in the inset to figure 8, all four resistivities collapse to essentially common 
form (that of [46] deviating just slightly at higher T). These differences are nonetheless 
potentially significant when comparing to theory, for which the magnetic contribution 
Pm&g(T) is required, obtained as in equation (2.1) by subtracting the resistivity of the 
non-magnetic LaCoIn^. The results in figure 8 show that the a factor in equation (2.1) 

- the relative weight of p(T) compared to that for the non-magnetic homologue - 
could vary by a factor of up to four or so, and is not therefore known with confidence. 

Figure 9 shows the experimental p^ g {T) (open circles) determined in [44] by 
subtracting the resistivity of LaCoIn§ from that for CeCoIn^; corresponding as such 
to a = 1 (or equivalently, if e.g. p(T) from [45] had been used instead, to an a of ~ 2). 
To compare to theory we consider the parameters e c = 0.5 and r\ = 0, with U = 3.75 
and V 2 = 0.8 (corresponding to an intermediate coupling strength U/irA ~ 1). The 
resultant theoretical p m &g{T) is compared to experiment in figure 9 (solid line), yielding 
a coherence scale of uo L = 60fC. It matches p^ g (T) for 15K < T < 100K, the deviation 
at low temperatures presumably reflecting the approach to the superconducting state. 
The deviation above T ~ 100/T also appears natural, since a direct determination of the 
CEF energy level scheme from inelastic neutron scattering [47] shows an excited level 
at 8.6 meV or ~ 100/T (with a second at a much higher energy, 24.4 meV), the extra 
conduction channel acting to diminish p^ g (T) more rapidly than the 1-channel theory. 

Our guess is that the latter inference is correct, at least qualitatively. A degree of 
caution is however required, since this is sensitive to a change in the value of a. To 
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Figure 8. Experimental resistivity oiCeCoIn§ measured by different groups: Ref [44] 
(circles), Ref [46] (squares), Ref [43] (triangles) and Ref [45] (crosses). Inset: showing 
collapse of experimental resistivities to common form on rescaling the y-axis alone. 
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Figure 9. CeCoIn^. Open circles denote the experimental p^i g {T) from [44], 
corresponding to a = 1 in equation (2.1) (dotted line shows the resistivity of LaCoIn^). 
The solid line shows the theoretical p ma g(T) for the bare parameters described in 
text, with coherence scale lol = 60K. Filled circles denote an experimental p^i g (T) 
obtained with a = 2.3 instead; the corresponding theoretical p m &g(T) (again with 
lul = 60K) is now shown as a point-dash line. Full discussion in text. 
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illustrate that, figure 9 also shows a new p^ g (T) (filled circles) obtained with a = 2.3 in 
equation (2.1). The theoretical p mSug (T) with the same coherence scale ujl = 60K, but 
with the overall y-axis naturally increased by factor of 2.3, now describes p^i g {T) very 
well for essentially all T > 15K. The quantitative influence of the extra conduction 
channel can thus be assessed with confidence only if the relative values of the resistivity 
for the magnetic and non-magnetic compounds are known accurately; although we add 
that the inferred low temperature coherence scale ujl — 60-K" is not itself sensitive to a. 

We turn now to the optical conductivity. Experimental results from reflectivity 
measurements [45] are shown in figure 10 (top panel), with the ^-dependence on a 
log scale up to w = 30cm -1 and a linear scale thereafter (separated by a vertical 
line). For uj < 30cm' 1 the results are extrapolated (dashed lines) towards the d.c. 
limit [45], as required for the Kramers-Kronig analysis of the reflectivity that leads to 
the experimental o~(uj; T) shown. The theoretical ct(uj; T) is shown in the lower panel, 
for the bare parameters specified above, with u L — 60K. 

The comparison between theory and experiment is at least qualitatively satisfactory, 
albeit not as good as for YbAl 3 or CeAl 3 (perhaps unsurprisingly given the quasi-2D 
nature of CeCoIn 5 ). The pseudogap at T = 10K is seen to lie at ~ 65cm -1 , and 
as in experiment shifts gradually to higher uj with increasing temperature and fills up 
progressively — on the scale of ~ 1 — 2oj l or ~ 60 — 120i^. The direct gap absorption 
in theory and experiment lies at Adh- ~ 600cm -1 ; and starts to lose spectral weight 
significantly for T > 50iC or so, i.e. on the scale of ujl itself. This suggests that 
CeCoIn*, is not in the strong coupling, Kondo lattice regime; since in strong coupling 
the direct gap/mIR absorption is significantly eroded for temperatures approaching the 
order of the direct gap itself (~ 600cm -1 or ~ 900K in the present case), see e.g. 
figure 12 of I. For CeCoIn 5 by contrast, significant thermal erosion is seen to occur 
for T > 50K ~ A dir /20, which behaviour is typical of intermediate coupling strengths. 
That also appears consistent with dHvA measurements [48], which yield a moderate 
effective mass m* in the range ~ 10 — 20. 

Two further points should be mentioned. First, the uj < 30cm -1 values of the 
theoretical c(uj; T) are clearly lower than the experimental extrapolations. As uj — > 0, 
the latter extrapolate to the d.c. conductivity obtained in [45]. The d.c. limit of the 
theoretical <r(0;T) by contrast gives the l/p mag (T) shown in figure 9 (solid line), and 
aside from the T = 292K case this agrees well with the experimental d.c. conductivity 
of [44] — which as seen from figure 8 differs by a factor of two or so from that of [45]. 
The issue here appears largely to be experimental, reflecting the significant difference 
between the resistivities of [44] and [45]. Second, the vertical arrows at uj ~ 250cm -1 
in the experimental o~(uj; T) of figure 10 indicate weak additional absorption that has 
been ascribed [45] to a Holstein band due to coupling to a bosonic mode. This is not 
of course included in the present theory, which thus shows somewhat less absorption in 
the region. 




Figure 10. Top panel: experimental optical conductivity of CeCoIn 5 at various 
temperatures [45]. Bottom panel: corresponding theoretical a(ui;T) obtained for the 
parameter set e c = 0.5, rj = 0, U = 3.75 and V 2 = 0.8, with uj l = 60K (~ 42cm" 1 ). In 
both cases the frequency axis is logarithmic up to w = 30cm _1 and linear thereafter 
(separated by a vertical line). Full discussion in text. 



Optical and transport properties of heavy fermions: theory and experiment 



20 



7. Conclusion 

We have here employed a local moment approach to the periodic Anderson lattice 
developed in I, to make direct comparison to d.c. transport and optical conductivities of 
CeB e , YbAlz, CeAl 3 and CeCoIn 5 . The Yb compound is a representative intermediate 
valence material, and the others typify heavy fermion behaviour, from the strongly 
correlated Kondo lattice regime appropriate to CeAl 3 and CeB 6 to what we believe 
is the somewhat weaker coupling case of CeCoIn^. In broad terms more or less all 
characteristic features of the optics and transport of these materials are captured; the 
natural exception, omitted from the model itself, being crystal field effects which may 
(or may not) show up in the experimental resistivity as a reduction below 1-channel 
behaviour at suitably high temperatures. The theory in general performs rather well 
quantitatively, and also captures notable features specific to individual systems — for 
example the existence of a low-frequency shoulder observed in the optics of YbAl 3 [31], 
or the absence of any significant direct gap/mIR absorption in CeAl 3 [35]. 

Minimalist though it is the underlying model, and theory for it, thus appear 
to provide quite a comprehensive and successful description of experiment. This we 
attribute in no small part both to the dominance of the local electron scattering inherent 
to the model itself, and the need to provide an adequate theoretical description of such 
on all experimentally relevant frequency and temperature scales. 
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